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Abstract 

About three million hadronic decays of the Z collected by ALEPH in the years 
1991 to 1994 are used to search for anomalous CV violation beyond the Standard 
Model in the decay Z — > bbg. The study is performed by analyzing angular 
correlations between the two quarks and the gluon in three-jet events and by 
measuring the differential two-jet rate. No signal of CV violation is found. For 
the combinations of anomalous CP violating couplings, h\ — \iAb9vb — hybgAb 

and hi = J hy b + h 2 Ab , limits of | hf, |< 0.59 and fc| < 3.02 are given at 95% 
CL. 
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1 Introduction 



The Standard Model predicts only negligible CP violating effects in decays of 
the Z into quarks or leptons |1. Therefore, a search for CP violation at the Z 
resonance is a test for physics beyond the Standard Model. Simple extensions 
H of the Standard Model with CP-odd and CP-even Higgs states are expected 
to contain CP violating terms. As couplings of Higgs particles are functions of 
the masses of the interacting particles, decays with heavy quarks deserve special 
interest. In this analysis |j] a search for CP violating anomalous couplings is 
performed, as proposed by the authors of jl], || . 



2 Theoretical Framework 

In [U CP-odd couplings are introduced in a model independent way using an 
effective Lagrangian. The Standard Model Lagrangian density is extended to 
include all CP-odd local operators that can be constructed with Standard Model 
fields, up to the mass dimension d = 6. Effects in Z — > bb from CP violating 
dipole form factors would require the measurement of the spin directions of the 
quarks. As for quarks no spin analyzers exist, the search for CP violation is 
restricted to the analysis of the CP-odd operator at the vertex Z — > bbg. In 
[0, fl it is shown that all CP-odd effects are proportional to the dimensionless 
coupling hi,: 

ht = h A b9vb - h V b9Ab (1) 

with 

- sva-dw cos$n/m^ 

"Ab/Vb = n,Ab/Vb 7 

where g s is the strong coupling and gvb, gAb are the Standard Model vector and 
axial vector couplings of the b quark to the Z. 

In a first study, CP-odd variables are analyzed, probing the coupling h\. To 
measure hb the CP-odd tensor Ty @] is used: 

\ \ K q X k q | J 

where k q , kq denote the normalized momentum vectors of the quark and anti- 
quark, respectively, and i, j are the cartesian coordinates with i, j = 3 defined 
to be along the beam axis. The tensor T'- is symmetric in i and j, traceless, 

and invariant when exchanging k q and kq. Therefore only the gluon jet has to 
be tagged. If CP is violated then (Ty) 7^ 0. The most sensitive observable is 
(T33) (Hi, it is related to h b j§ by: 

(T> 3 ) = h b Y> . (3) 
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The sensitivity Y' is a constant derived by integration of the experimentally 
accessible phase space. 

In a subsequent analysis, additional contributions from anomalous couplings 
to the three-jet rate in the process Z — ► bbg have been searched for. Such 
additional contributions are proportional to the combination h^, with 

K = y/h 2 Vb + h 2 Ab , (4) 

and would manifest themselves in a higher value of the strong coupling constant 
for b quarks, a b s (M^) (for a complementary but related study see Q). The 
analysis of is based on the differential two-jet rate D2, which is the normalized 
distribution of the event shape variable Y 3 , D^iY^) = 1/a da(Y 3 )/dY 3 , with Y 3 
being the y cu t at which an event changes its classification from three jets to 
two jets. Additional couplings at the vertex Z — > bbg could also contribute to 
the partial width T^, thus possibly enhancing R b beyond its Standard Model 
value g|. By attributing the observed deviation of R e b xp = 0.2209 ± 0.0021 
from the Standard Model expectation i?/, = 0.2155 ± 0.0005 || entirely to these 
anomalous couplings, a value of ht — 1.93 ± 0.75 can be calculated. 



3 Data Reduction 

The ALEPH detector and its performance are described in detail in jj. The 
selection of hadronic events is detailed in || . Jets are defined using the JADE 
or Durham algorithms ||| . In the h b analysis the JADE algorithm with a fixed 
cut-off value of y cu t = 0.03 is used. In the hi analysis both jet schemes are 
employed to measure the differential two-jet rate D2(Y 3 ). The b events are 
selected exploiting the sizeable impact parameters due to b decays. The b tagging 
algorithm is described in detail in reference || . The probability for each charged 
particle to originate from the primary vertex is calculated using the track impact 
parameter. They are combined into probabilities V associated to sets of tracks 
like jets (Vj) or whole events (Vg). In both analyses b events are selected with 
a cut on the event probability Ve, resulting in b purities of more than 85%. 

4 Analysis of Anomalous Couplings from Z — > 
bbg Topologies 

In the study of CP-odd contributions three-jet events are selected as described 
in ||. The energies of the three reconstructed jets are ordered, j > £j e t 2 > 
-Ejg^g. In the b sample the gluon jet candidate is selected using the probability 
Vjoi the individual jets: from the two lower energy jets, the jet with the higher 
Vj is chosen as the gluon jet candidate. 



2 



The purity of the sample is estimated using simulated events generated with 



JETSET 1 10 . The matching procedure of partons to jets found after recon- 
struction is detailed in Sj. A purity of 73.6% ±0.4 stat ±0.9 sys t is achieved. The 
efficiency, defined as selected bbg events with successfully tagged gluon com- 
pared to the total rate of Z — > bbg events, is about 19%. The measurement is 
performed on the data collected until 1994. After applying the data selection, a 
sample of 85342 bbg candidates remains for the Tg 3 measurement. The following 
mean value of (T^ 3 ) is obtained: 

<^ 3 ) = (-0.5±3.7 stat )xl(r 3 

The measurement is also performed on a sample enriched with light quarks 
to ensure that no effect due to the selection mechanism is present and to check 
the assumption that no CV violating effect exists in events with light quarks. 
Light flavours are selected with the cut Vs > '-V- The gluon jet in this sample 
is defined to be the one having the lowest energy in the event. The following 
mean value of (T33} is obtained: 

(^ 3 ) = (-0.9±3.1 stat )xl0- 3 

Obviously no significant discrepancy between the two completely independent 
samples is found. Hence it is concluded that no significant fake CV violating 
effect due to the selection is present. 



4.1 Systematic Errors 

To estimate the effect of the selection cuts, all cuts are varied in a wide range. 
Only those applied to select well-defined three-jet configurations (y C uti a Pl a_ 
narity, jet energies, multiplicities and angles) have a sizeable influence on the 
selected sample. Therefore, the systematic error is estimated using these cuts 
only j|. The CV invariance of the fotag is checked by dividing the three-jet 
sample into disjoint subsamples with different Vs values. Using these subsam- 
ples, a sample of (T33} values is measured to test the bias of the 'Pg cut. The 
(T33) values are found to be well compatible with statistical fluctuations. The 
uncertainty is estimated from a Gaussian fit to the distribution. The CV sym- 
metry of the detector is estimated by measuring (T33) on a sample of track pairs 
in hadronic events || . The resulting asymmetry is compatible with zero within 
one standard deviation. The quality of the reconstruction of the jet direction is 
estimated by measuring the difference of the (T33) at the parton level and after 
the full detector simulation. The error is derived from the width of the differ- 
ence distribution. A more detailed description of the systematic uncertainties 
is given in H . 
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4.2 Determination of hb 

To extract the coupling hb from the measurement of (T33) the effective sensi- 
tivity Y' has to be calculated. This has been done by means of a Monte Carlo 
generator fuJI , which includes the CP violating couplings. The calculated sen- 
sitivity Y' = —0.0167 ± 0.0002 s ^ a ^ is constant, showing no dependence on h\ 
[gj . The systematic errors on the sensitivity stem from the errors on the tagging 
purity and from the influence of the b quark mass. The largest error is caused by 
the dependence of the sensitivity on the b quark mass |T^| . Taking into account 
the total systematic uncertainty || the sensitivity is given by 

Y' = -0.0167 ± 0.0002 stat ± 0.0015 syst . 

4.3 Results of the h b Analysis 

Taking the systematic errors into account the measurement of the CP-odd ob- 
servable in a sample of Z — > bbg events yields: 

(T 3 ' 3 ) = (-0.5 ± 3.7 stat ± 3.3 syst ) x 10~ 3 . 

The measurement is consistent with (T33) = 0. The size of the CP-odd coupling 
hb is extracted using eq. |[ 

h b = 0.03 ± 0.22 stat ± 0.20 syst . 

From this measurement a limit on the coupling of | hb \< 0.59 (95% CL) is 
derived. 

5 Measurement of Anomalous Couplings from 
the Differential Two-Jet Rate 

Of all the event shape variables studied in || , only for Y3 theoretical calculations 
including anomalous couplings are available jl2| . Therefore, the second analysis 
concentrates on possible additional contributions to the differential two-jet rate 
-^2(^3) due to the anomalous coupling h^. The theoretical prediction for D2(Y 3 ) 
is given by 

D h 2 (Y 3 ) = K^) + hfC'(Y 3 ) 

Z7T \ 

Ai(Y 3 )2nb Q \n(^j+B'(Y 3 ) 

where 60 = (33 — 2n/)/127T, \x is the renormalization scale, and the coefficients^' 
and B' have been computed to second order of perturbative QCD pa]. Here 
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ft| 3> /i* dsc is assumed. Hence the prediction for udsc quarks is entirely fixed by 
the Standard Model and does not depend on h* . In contrast, D2 for b quarks 
receives contributions from new physics and the prediction is modified by an 
additional term proportional to Note that in the coefficients A' and B' the 
total cross section needed for normalization is changed if anomalous couplings 
are present. The coefficient C has been calculated to leading order |L2[ , 
A theoretical description i?^h Q f ^ ne observable can be derived: 

Rth p b p b 2 bi&g {Y z ) + (1 - n)^ dscMag (y 3 ) 

R b D b 2 mc \Y 3 ) + (1 - R b )Df sc ' mcl (Y 3 ) 

where P b is the purity of the lifetime-enriched b sample, R b is the fraction of Z's 
decaying into b quarks and D\ tag stands for the distribution of a flavour q in a 
sample of type tag. These distributions are constructed from the parton level 
predictions and have to be corrected for the following effects: mass corrections 
for b quarks, initial and final state radiation, hadronization effects, the detector 
acceptance, the influence of the detector resolution, and the tagging bias. After 
applying these corrections to i? , a value of hi is extracted from a binned 
least-square fit to the data. 



5.1 and Systematic Errors 

The fit function i?^ n depends explicitly on h%, a s and /i. The value of the 
strong coupling constant is set in the fit function to a s (M|) = 0.II8 and is 
varied by ±0.007. The renormalization scale is set to 15 GeV. This symmetrizes 
the error from the scale uncertainty, which is varied from \i = m b to fi = Mz- 
The fit range is chosen to optimize the sensitivity and to guarantee a good 
perturbative description. To take into account systematic uncertainties due to 
the limited fit range the fit is repeated with a range modified by two bins. 
Jets are reconstructed using charged tracks and neutral calorimeter objects. 
The uncertainty of the reconstruction procedure is estimated by repeating the 
analysis with charged tracks only. The b sample is enriched by means of a 
lifetime tag, which leads to a distortion of the differential two-jet rate of less 
than 10%. Therefore, a correction is elaborated using full detector simulation. 
The stability of these corrections is checked by varying the lifetime cuts in 
the data resulting in a change of the purity of the sample of 10%. The same 
cuts are applied to the simulated data and the corrections for the tagging bias 
recalculated. The b quark fragmentation is described by the fragmentation 
function of Peterson et al. |l^| . The main parameter of this function is e b , 
measured to be e b = (3.2 ± 1.7) x 10~ 3 fl^j . Monte Carlo simulations with a 
corresponding range of values are done to study the effect on hi . Another source 
of uncertainty is related to mass corrections. These have been calculated in [[l6f 
at tree level. These calculations are only complete to 0(aj) and are applied to 
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the coefficient A'. The uncertainty on the b quark mass is set to 0.5 GeV/c 2 
and the corresponding correction recalculated. In order to account for missing 
higher orders, the available 0(aj) four-jet computation is used for correction as 
well and the difference to the 0{ar) result is taken as systematic error. Finally, 
the parameters Rt and Pb are varied in the fit function within their errors and 
an error on the normalization is derived. The typical systematic uncertainty of 
the cluster algorithms is Ahl = 0.38 ||. 

5.2 Hadronization Model Uncertainty and Results of the 
Analysis 

The main theoretical uncertainty is caused by the hadronization. The results 
achieved with four different generators are compared in table |l|: the matrix 
element with string fragmentation (ME), the parton shower (Qq — 1 GeV, being 
the cut-off of the parton shower) as implemented in JETSET [|l0| with string 
fragmentation (PS), the model of cluster fragmentation in HERWIG |l7| (HW) 
and the dipole cascade model implemented in ARIADNE jl8) (AR) . Taking into 
account the correlations between the two algorithms, a combined result for each 
hadronization model can be derived which minimizes the total statistical error. 



Model 


h* 


ME 


1.05 ±0.26 ±0.40 


AR 


0.79 ±0.30 ±0.40 


PS 


1.75 ±0.16 ±0.35 


HW 


1.79 ±0.16 ±0.36 



Table 1: Results on hi from the least-square fit for four different hadronization 
models. The first error is the statistical one, the second the total systematic 
uncertainty. 



After averaging the combined results of the different hadronization models 
the following result ist obtained: 

ht = 1.34 ± 0.22 stat ± 0.38 syst ± 0.50 hadr . 

where the first error is the statistical one, the second error contains all sys- 
tematic uncertainties but the hadronization model and the last one reflects the 
systematic uncertainty due to the hadronization model itself. The upper limit 
is derived by adding linearly the systematic and statistical errors: 

h* h < 3.02 (95% CL) . 
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6 Conclusions 



A search for CV violation beyond the standard model in the decay Z — » bbg 
has been performed. Two combinations of CP-odd couplings, namely hi and 
hb have been analyzed. No evidence for CP-odd couplings is found in both 
analyses. The derived limit on is consistent with the value calculated from 
the i?b measurement. 

Using eq. [l] and eq. ^, the two measurements presented in this paper can be 
used to constrain the couplings h,Ab, hvb- I n n g- the 95% CL limits of both 
measurements on fiAb and hyb are shown, being well consistent with each other. 




-4-2 2 



Figure 1: Combined results. The shaded areas depict the constraints of the 
measurements on the couplings h Ab and hyb- 
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